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NATIONAL .ADVISORY. .COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1406 


ON THE USE OF THE HARMONIC LINEARIZATION METHOD IN THE 
AUTOMATIC CONTROL THEORY* 

By E. P. Popov 

The method of harmonic linearization (harmonic "balance) , first 
proposed by N. M. Krylov and N. N. Bogolyubov (ref. l) for the approxi- 
mate investigation of nonlinear vibrations, has been developed and re- 
eived wide practical application to problems in the theory of automatic 
control (refs. 3 to 6) . Recently, some doubt has been expressed on the 
legitimacy of application of the method to these problems, and assertions 
were made on the absence in them of a small parameter of any kind. Never- 
theless, the method gives practical, acceptable results and is a simple 
and powerful means in engineering computatations . Hence, the importance 
of questions arises as to its justification. The underlying principle of 
the method is the replacement of the given nonlinear equation by a linear 
equation. In establishing the method, a small parameter is considered 
whose presence makes it possible to speak, with some degree of approxi- 
mation, of the solution of this new equation to the solution of the given 
nonlinear equation. In an article by the author (ref. 7), certain con- 
siderations were given on the presence of the small parameter, but this 
question has not as yet received a final answer. In the present report, 
a somewhat different approach to the problem is applied that permits: 

(a) establishing, in the clearest manner, the form of the presence of 
the small parameter in nonlinear problems of control theory, solvable by 
the method of harmonic linearization; (b) connecting it with previous 
intuitive physical concepts (with the "filter property") and extending 
the class of problems possessing this property; and (c) discussing various 
generalizations of the method. 

The free motion (transition process and autovibration) for a very 
wide class of nonlinear systems of automatic control (ref. 7) are 
described by differential equations of the form 

Q(p)x + R(p)F(x,px) =0 (p = (1) 
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where Q(p) and R(p) are polynomials of "any degree, of which some re- 
quired properties will he established in the following paragraphs, and 
F(x,px) is a given nonlinear function possible only with respect to as- 
sumptions of the most general character. However, in problems of the 
theory of control, no assumptions must be made as to the smallness of 
the nonlinear function F(x,px) or to its small difference from a linear 
function. In order to render explicit the form in which it would be 
possible to write a small parameter in equation (l), we shall proceed as 
follows . 

Let equation (l) have a periodic solution or a solution approximately 
periodic differing slightly from the sinusoidal. We write this solution 
in the form 


x = x* + ey(t) x* = a sin cot (2) 

where e denotes a small parameter, and y(t) denotes an unknown bounded 
function of time. In the case of the existence of a periodic solution, 
we write 


ey/t) = \ sin(koot + q> k ) (3) 


We represent the given nonlinear function F(x,px) in the form 

F(x,p) = F(x*,px*) + [F(x* + ey, px* + epy) -F(x*,px*)] (4) 

Expanding the two components separately in a Fourier series, we obtain 


F(x,px) 


*30 


(f + ip) 


sin cot 


+ ZZ F k + 

k=2 


* 12 % 

k=0 


( 5 ) 


where q q. A, and B are the coefficients of the initial term, the sine, 
and cosine terms, rexpectively , the cosine being replaced by — p sin cot 

of the expansion of the function F(x*,px*) in a Fourier series. ZF k 
denotes all the higher harmonics of the expansion of F(x*,px*) in a 
Fourier series (they must not be considered small, since the nonlinearity 
is not small) where we write 

F k = fe k sin ( la0fc + ^k) (k = 2, 3, . . ., ®) (6) 


eZ^^. denotes all the terms of the expansion in a Fourier series of all 
the expressions shown in brackets in formula (4) . This entire expression 
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is written with, a small parameter which, according to equation (4) , is 
small if the derivatives Sf/Sx and dF/5px are finite. This expression 
is also computed as small in the case of certain discontinuous nonlinear 
characteristics (e.g., the Raleigh type where the preceding derivatives 
at the points of discontinuity are delta functions) . We may write 


e<i> = 6c. sin(kcut + ) (k = 0, 1, 2, . . . ,*) 


(?) 


We substitute equations (2) and equation (4) in the given equation (l) , 
so that 

Q(p)x* + Q(p)ey + 

V oo OP 

R(p)q n + R( P )(a + | p) sin cot + R(p)e^ F k + = 0 (8) 

\ ' k=? k=0 


Since the equation must be satisfied identically, we separately 
equate to zero all coefficients with the same order harmonic. We note 
that formula (8), in the case of the existence of a periodic solution, 
is exact. 

From the equating of the zeroth harmonics of equation (8), there is 
obtained with an accuracy up to e the relation 


'2jt 


-I F(a sin u, acu cos u) du = 0 (u = cut) 


( 9 ) 


which is a certain general requirement for F(x,px). 

From the equating of the first harmoics of equation (8) , taking 
account of equations (2) and equation (7), we have 


a sin cut 


= - V 


+ B* 




+ T + P) - 


6 C-i 


R(icu) 

Q(icu) 


sin (cut + ^ 2 . + P) 


( 10 ) 


where x aud (3 are arguments of the expressions A + iB and 
R ( icu) / Q( icu) . On the basis of the exact equation (10), we obtain the 
following approximation (with an accuracy up to e) : 


a 




A c + B 


R 

(icu) 

Q 

[icu) 


x + 3 * it 


(ID 
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It is here assumed that the polynomial Q(p) in equation (l) does not 
have purely imaginary roots. 

From the equating of the higher harmonics of equation (8), consider- 
ing equations (3), (6), and (7), we have 

eajj sin(kmt + cp^.) 




||^|sin(kmt + t k + P k )- 6c k||[iS]'| sin ( tot + + Pk) ( 12 ) 


where 0 k denotes the argument of the expression R ( ikn>) / Q( ikn>) . 


It is thus seen that if b^. is not small, the magnitude 
| R(ikcu)/Q(ikm)| should he of the order of e. The last component in 
equation (12) will then he of the order From the exact equation 

(12), we obtain the following approximate equation (with an accuracy up 
to e) : 


ea k 


R(ikco) 
Q( ikes) 


Comparing formulas (13) and (ll) 
to have in the solution (see eqs. (2) 



*k + ?k ” * 

we see that, for example, 
and (3)) 


(13) 
the wish 

(14) 


leads to the need of satisfying, in the given equation (l), the following 
requirement : 



R(ikcio) 
Q( ikoo) 


<< 


(A 2 + B 2 ) 


R(ico) 


Q( ico) 


(15) 


where its satisfying in the concrete system can he checked after co is 
obtained. The degree of Q(p) should, in any case, he higher than that 
of R(p) . A particular case of the general expression (15) is intuitively 
the earlier introduced "filter property" of the linear part. 

Thus, condition (15) has been obtained and must he satisfied by the 
coefficients of the given differential equation (l) in order that a pe- 
riodic solution, if it exists, may he approximately determined in the 
form of sinusoids in the presence of a "strong" nonlinearity of F(x,px) . 
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The equation Tor its approximate determination according to equation (a) 
with the substitution of sin cot = x*/a, assumes the form 


[q(p) + R(p) (q + 7J7 p)] x* = o 


(16) 


where 




aco cos u)sin u du 




sin u, aco cos u)cos u du 


(17) 


The replacement of equation (l) by equation (16) with its subsequent 
investigation by the linear methods is called harmonic linearization. 

This is formally equivalent to following the mode of writing the initial 
equation (l) with a small parameter, so that 

Q(p)x + R(p) (q + ^ p) x + £ f(x,p) = 0 (18) 

where, according to equation (8), we may write 


ef(x,p) = R(p) 


£ 

k=2 


F k +s 


-(q + ^ p) £ y 


k=0 


where the first term of the three terms shown in brackets, taken sepa- 
rately, is not small. The smallness of the function ef(x,p) is acquired 
only with the factor R(p) owing to the property of equation (15) . 

There has thus been found the form of the presence of the small 
parameter s in the equations of nonlinear systems of automatic control 
required for the application of the method of harmonic linearization 
(this also refers to the first approximation of the method of small pa- 
rameter which, according to reference 7, coincides with the given method) . 
The writing of equation (l) with small parameter in the form of equation 
(18) is valid in the region of the existence of equation (2). 

For definiteness we shall assume that the polynomials R(p) and 
Q(p) are such that the characteristic equation (eq. (16)) has all positive 
coefficients, and that the quotient of the division of the entire left 
side of equation (16) by p 2 + cn 2 satisfies the criterion of Hurwitz. 

It may then be said that the given nonlinear system (eq. (l)) for the 
existence of the periodic solution (eq. (2)) is close to the equivalent 
linear system which is on the boundary of stability, except for the 
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circumstance that the linear system is different for various periodic 
solutions (because of a change in the coefficients q and q^ from one 
periodic solution to another, since q and q^ depend on the amplitude 
a and the frequency ao of the required solution). Thus, evidently, the 
equation of the first approximation (eq. (16)) is nonlinear if one speaks 
of the combination of all possible periodic solutions for different values 
of the coefficients of the polynomials R and Q (i.e., for different 
parameters of the system) , and is linear only for each given periodic 
solution (along the given solution) . Such is the property of the method 
of harmonic linearization in its application to systems with "strong" 
nonlinearity . 

By writing the given equation (l) in the form of equation (8) all 
the variants of the method of harmonic balance in control theory and all 
generalizations of the method also become easily understandable. 

A similar reasoning may also be developed for the application of the 
method of -N. N. Bogolyubov (ref. 2) to the investigation of transient 
vibrational processes in certain nonlinear systems of control. 
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